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SUMMARY 


The  small  separation  bubbles  which  form  near  the  leading  edge 
of  airfoils  prior  to  the  onset  of  leading-edge  stall  have  been 
analyzed  in  detail,  including  the  effects  of  viscous-inviscid 
interaction.  The  separated  laminar  shear  layer,  transitional 
flow  and  turbulent  reattaching  flow  are  represented  by  an 
integral  formulation,  h correlation  of  local  shear-layer 
parameters  has  been  developed  for  determining  the  onset  of 
transition  in  the  laminar  shear  layer.  Solutions  are  obtained 
using  an  iterative  procedure,  with  strong  interaction  effects 
limited  to  the  immediate  vicinity  of  the  separation  bubble. 
Results  obtained  for  specific  airfoils  are  in  good  agreement 
with  wind  tunnel  measurements.  The  method  was  used  to  investi- 
gate the  mechanism  for  bubble  bursting.  Results  indicate  that 
reseparation  of  the  turbulent  boundary  layer  downstream  of 
reattachment,  rather  than  failure  of  the  shear  layer  to  re- 
attach, causes  bubble  breakdown. 
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1.  INTRODUCTION 


The  regions  of  separated  flow  which  form  on  airfoils  govern 
the  airfoil  stall  characteristics.  The  nature  and  extent  of 
these  regions  are  determined  primarily  by  the  airfoil  shape, 
incidence  and  Reynolds  number.  Of  concern  here  are  the  small 
separation  bubbles  which  form  near  the  leading  edge  of  air- 
foils of  moderate  thickness  ratio  (.09  to  .15)  at  chordal 
Reynolds  numbers  in  the  range  from  one  to  ten  million.  The 
occurrence  of  what  is  termed  leading-edge  stall,  characterized 
by  an  abrupt  loss  in  lift  and  increase  in  drag,  can  be  attri- 
buted to  the  sudden  breakdown,  or  bursting,  of  the  leading  edge 
bubble  (Reference  1)  . 

The  flow  in  the  vicinity  of  the  leading  edge  of  an  airfoil 
subject  to  leading-edge  stall  is  as  sketched  in  Figure  1.  The 
laminar  boundary  layer,  extending  from  the  stagnation  point 
over  the  leading  edge,  separates  just  downstream  of  the  point  of 
minimum  pressure.  Transition  to  turbulent  flow  occurs  in  the 
free  shear  layer  a short  distance  downstream  of  the  separation 
point.  The  flow  then  reattaches  to  the  airfoil  surface,  with 
a turbulent  boundary  layer  extending  from  the  reattachment 
point  to  the  trailing  edge.  If  the  angle  of  attack  of  the 
airfoil  is  increased,  the  bubble  moves  closer  to  the  leading 
edge  and  becomes  slightly  shorter.  The  bubble  has  almost  no 
effect  on  integrated  loads,  because  it  is  never  more  than  a 
few  percent  of  the  chord  in  length.  In  the  immediate  vicinity 
of  the  bubble,  though,  there  is  strong  interaction  between 
the  viscous  and  inviscid  flows. 

The  specific  mechanism  for  bubble  breakdown  is  not  presently 
known.  It  has  been  postulated,  though,  that  there  is  a physical 
limitation  in  the  amount  of  pressure  recovery  possible  in  the 
turbulent  shear  layer,  so  the  bubble  bursts  when  the  limit 
is  exceeded  and  the  shear  layer  fails  to  reattach.  Alternatively 
it  has  been  suggested  that  stall  results  from  separation  of  the 
turbulent  boundary  layer  just  downstream  of  reattachment  (Refer- 
ence 2).  This  study  was  undertaken  to  provide  a tool  for  invest! 
gating  the  specific  mechanism  for  breakdown  of  the  leading-edge 
bubble  and,  ultimately,  for  accurately  predicting  the  onset  of 
leading-edge  stall. 

Leading-edge  bubbles  have  been  the  subject  of  numerous  studies. 
Reviews  of  this  work  are  given  in  References  2 and  3.  The 
primary  difficulties  in  treating  the  problem  analytically  derive 
from  the  interaction  between  the  viscous  and  inviscid  flows  and 
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the  coupling  between  the  interaction  and  transition  from 
laminar  to  turbulent  flow  in  the  free  shear  layer.  Analyses  have 
been  carried  out  using  semi-empirical  formulations  which  do 
provide  fairly  good  qualitative  agreement  with  tests  in  pre- 
dicting stall  onset  (References  4 and  5) . However,  the  modeling 
has  not  been  adequate  for  analyzing  the  details  of  the  flow 
in  and  near  the  bubble. 

More  recently,  separation  bubbles  which  occur  near  midchord 
on  relatively  thick  airfoils  at  zero  incidence  were  analyzed 
numerically  using  a finite-difference  method  (Reference  6). 

While  these  bubbles  are  about  ten  times  longer  than  leading- 
edge  bubbles,  their  structure  is  quite  similar,  with  transi- 
tion occurring  in  the  free  shear  layer.  The  results  of  that 
analysis  are  in  excellent  agreement  with  the  flow  measure- 
ments reported  in  Reference  7.  It  should  be  noted,  in  particular, 
that  in  Reference  6 interaction  was  assumed  to  be  limited  to 
the  vicinity  of  the  bubble.  Also,  the  validity  of  the  boundary- 
layer  approximation  for  analyzing  the  bubble  was  verified  by 
direct  comparison  with  a solution  using  the  complete  Navier- 
Stokes  equations. 

In  this  study,  the  separated  and  reattaching  shear  layer  in  a 
leading-edge  bubble  were  analyzed  using  an  integral  formula- 
tion, assuming  the  boundary- layer  approximation  is  applicable. 
Interaction  between  the  viscous  and  invescid  flows  in  the 
vicinity  of  the  bubble  was  taken  into  account  through  an 
iterative  procedure.  The  method  of  analysis  is  outlined  briefly 
in  the  next  section.  Derivation  of  the  formulations  employed 
and  results  of  analyses  for  different  airfoil  sections,  angles 
of  attack  and  Reynolds  numbers  are  presented  in  Reference  8. 

The  method  developed  for  analyzing  leading  edge  bubbles  was 
employed  to  investigate  the  specific  mechanism  for  bubble  burst. 
Details  of  that  analysis  are  given  in  Section  3. 
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2.  METHOD  OF  ANALYSIS  OF  LEADING-EDGE  BUBBLES 


Inviscid  Flow  Representation 

Assv.uiing  the  flow  to  be  two-dimensional,  let  Uo(x)  denote  the 
magnitude  of  fluid  velocity  at  the  airfoil  surface  that  would 
result  in  the  absence  of  viscous  effects,  where  x is  a co- 
ordinate measured  along  the  airfoil  surface.  The  flow 
component  tangent  to  the  airfoil  surface  at  the  interface 
of  the  viscous  and  inviscid  flows  is  written  in  the  form 

Ue(x)  = Uq(x)  1 + u (x)J  (1) 

A previously  developed  digital  computer  program  is  employed 
in  the  analysis  to  direct  the  computation  of  Uq,  given  the 
airfoil  shape  and  angle  of  attack,  using  a source  distribu- 
tion on  the  airfoil  surface  and  a vortex  distribution  on  the 
chord  line. 

The  analyses  of  the  boundary  layer  and  shear  layer  provide  the 
flow  inclination  Vq/uq  at  the  interface  of  the  viscous  and 
inviscid  flows,  which  can  be  related  to  u as  follows.  Inter- 
action is  taken  to  occur  on  the  interval  x;^  ^ x ^Xg,  over 
which  it  is  assumed  the  surface  curvature  is  negligible.  The 
perturbation  to  the  inviscid  flow  is  derived  from  a potential; 
the  potential  is  formulated  from  a source  distribution  on  the 
airfoil  surface. 


With  m(x)  = Ve^e  Ve/u^, 

and 

assuming 

u < < 1, 

it  is  found  that 

1 

u(x)  - 

Tf 

f 

m(  f 

) d? 

(2) 

T 

^A 

X - 

f 

Viscous  Flow  Representation 

The  flow  in  the  free  shear  layer  and  the  boundary  layer  in 
the  vicinity  of  the  separation  bubble  are  represented  using 
the  integral  formulation  developed  in  Reference  9 for  analyzing 
supersonic  separated  and  reattaching  laminar  flows  involving 
strong  interaction  with  the  inviscid  flow.  The  relations 
have  been  generalized  to  account  for  continuous  transition 
from  laminar  to  turbulent  flow  in  the  free  shear  layer.  Both 
the  momentum  integral  (zeroth  moment)  and  first  moment  of 
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momentum  of  the  boundary  layer  equations  are  used,  so  that  for 
laminar  or  fully  turbulent  flow  the  velocity  profiles  can  be 
characterized  by  a single  parameter  which  is  not  related  to 
the  local  pressure  gradient.  The  family  of  similar  solutions 
for  reversed  flow  found  by  Stewartson  (Reference  10)  is 
employed  for  analyzing  the  free  shear  layer.  Turbulence 
production  is  introduced  using  an  exponentially  increasing 
intermittancy  function,  with  the  constant  in  the  exponent 
determined  from  measurements  taken  in  a free  shear  layer  under- 
going transition  (Reference  11) . 

A coupled  pair  of  first-order,  ordinary  nonlinear  differential 
equations  was  derived  from  the  momentum  equations.  They  have 
been  formulated  in  such  a way  that  they  can  be  integrated 
continously  downstream,  starting  in  the  laminar  boundary  layer 
and  continuing  through  the  separation  point,  transition  in  the 
laminar  shear  layer  and  reattachment  of  the  turbulent  shear  layer. 
The  dependent  variables  are  the  displacement  thickness  §*  and 
a parameter,  denoted  a,  characterizing  the  velocity  profile. 

That  parameter  is  defined  as  follows: 


a 


, attached  flow; 

0 


(3a) 


separated  flow. 


(3b) 


The  integral  across  the  boundary  layer  of  the  momentum  equation 
and  of  that  equation  multiplied  by  u then  give,  under  the  afore- 
mentioned assumptions,  a pair  of  differential  equations  of  the 
form 


d 5 * 
dx 

(4) 

da 

dx 

= (fa  + ga  Ue  ) / Ue  + h^ 

(5) 

where  the  f's,  g's  and  h's  are  nonlinear  functions  of  S*  and 
a.  With  Ue(x)  specified,  Eqs . (4)  and  (5)  can  be  integrated 

numerically  on  x to  obtain  S'*  and  a. 

To  complete  the  basic  formulation  of  the  viscous  flow,  the 
continuity  equation  provides  the  following  relation  for  flow 
inclination  at  the  edge  of  the  layer; 


( 


!i 


IJ 
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A relation  of  the  form 


was  postulated  for  correlating  transition  onset  location  in 
the  bubble,  where  subscripts  s and  t refer  to  separation  and 
transition  onset,  respectively.  Using  the  measured  pressure 
distributions  given  in  Reference  7,  Eqs . (4)  and  (5)  were 
integrated  through  separation  and  into  the  laminar  free  shear 
layer  to  provide  streamwise  variations  of  a ^ . A total  of 
15  cases  were  analyzed,  including  leading-edge  bubbles  on 
two  different  airfoils,  each  at  two  different  angles  of 
attack,  for  free-stream  Reynolds  numbers  ranging  from  1.5  x 10^ 
to  6 X 10^.  For  each  case,  the  measured  location  of  transition 
given  in  Reference  7 was  used  to  obtain  the  appropriate  value 
of  a^  St/S%  • The  results  obtained  are  plotted  in  Figure  2. 
While  there  is  some  scatter,  particularly  at  the  higher  Reynolds 
numbers,  the  data  are  still  well  correlated  by  the  postulated 
relation.  The  line  drawn  through  the  points  is  a plot  of  the 
simple  relation 


which  is  seen  to  provide  a quite  good  approximation  to  the 
derived  correlation.  This  equation  was  used  in  the  bubble 
analyses  to  locate  transition  onset. 

Iteration  Procedure 


An  iteration  step  is  begun  by  integrating  Eqs.  (4)  and  (5)  to 
obtain  S*  and  a as  a function  of  x,  given  Ue  and  u4*  The 
step  is  completed  by  obtaining  revised  estimates  of  Ue  and  Ue 
as  dictated  by  the  relations  governing  the  viscous-inviscid 
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CORRELATION  OF  SHEAR  LAYER  PARAMETERS 
AT  TRANSITION  ONSET  IN  A LEADING-EDGE  BUBBLE 


FIGURE  2 
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interaction.  The  obvious  procedure  of  simply  substituting 
the  variation  of  Vg/ue  obtained  from  the  viscous-flow  analysis, 
Eq.  (6),  into  the  integrand  of  Eg.  (2),  was  found  to  be 
unsuitable.  The  initial  estimate  for  the  Ug  variation  inevitably 
produces  small  but  physically  unrealistic  excursions  in 
in  the  immediate  vicinity  of  transition  which  cause  rapid 
divergence  of  results  for  successive  iterations.  The  following 
procedure  was  therefore  devised. 

The  interaction  is  introduced  using  the  differential  equation 
for  Eq.  (4),  as  a link  to  the  variation  of  obtained 

in  the  analysis  of  the  viscous  flow.  Specifically,  a and  f*, 
and  hence  the  coefficients  f fa»  etc.,  are  regarded  as 

known  in  Eq.  (4),  their  variation  having  been  determined  from 
the  viscous-flow  analysis,  while  d S */dx,  dug/dx,  and  Ug  are 
regarded  as  unknown.  By  combining  Eqs . (1),  (2)  and  (6)  in 

Eq.  (4),  one  then  obtains  a linear  integro-dif ferential  equation 
for  the  flow  inclination  m(x) . This  equation  is  solved  using 
techniques  analogous  to  those  employed  in  thin-air foil  theory. 

A weighted  average  of  the  previous  and  derived  variations  of  Ug 
is  then  employed  to  begin  the  next  iteration  step. 

Representative  Results 

Results  of  a bubble  analysis  for  the  case  of  a modified  NACA 
0010  airfoil  at  8 degrees  angle  of  attack  and  chordal  Reynolds 
number  of  2 x 10^  are  shown  in  Figure  3,  where  variations 
between  separation  and  reattachment  of  Ug,  S , S*  and  the 
ordinate  where  u = 0 are  plotted.  The  agreement  between 
computed  and  measured  variations  of  u^  and  the  locations  of 
reattachment  (measured  location  is  marked  by  a small  arrow 
on  the  abscissa)  are  seen  to  be  very  good.  Further  results 
are  given  in  Reference  8. 
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3.  ANALYSIS  OF  BUBBLE  BURST  MECHANISM 


An  analysis  was  undertaken,  using  the  method  outlined  in  the 
previous  section  as  the  basic  tool,  with  the  aim  of  determining 
whether  bubble  burst  can  be  attributed  to  either  failure  of  the 
shear  layer  to  reattach  or  to  reseparation  downstream  of  re- 
attachment. The  NACA  0012  airfoil  section,  which  is  generally 
accepted  as  being  subject  to  leading-edge  stall,  and  for  which 
there  is  a large  body  of  data  concerning  stall  available,  was 
selected  as  the  specific  subject  for  analysis. 

As  a first  step,  the  potential  flow  and  loading  for  incipient 
stall  was  defined.  Data  on  maximum  lift  coefficient  max 
was  obtained  from  References  12  through  15,  and  the  plot  of 
C/jjiax  • chordal  Reynolds  number  shown  in  Figure  4 was  generated. 
A Reynolds  number  of  2 x 10^  was  selected  for  detailed  analysis, 
as  that  value  is  well  within  the  Reynolds  number  range  over 
which  leading-edge  stall  takes  place  (roughly  10^  to  6 x 10^  - 
see  Reference  2).  From  Figure  4,  it  is  seen  that  at  this 
Reynolds  number  bubble  burst  should  occur  at  a lift  coefficient 
of  1.37. 

Addressing  first  the  question  of  whether  burst  is  due  to  failure 
to  reattach,  analyses  of  leading-edge  bubbles  on  a 0012  airfoil 
at  a Reynolds  nurrOser  of  2 x 10^  were  carried  out  using  the 
procedures  described  in  the  previous  section,  for  lift  coefficients 
of  1.27  (i.e.,  somewhat  below  stall),  1.37  (incipient  stall), 
and  1.57  (well  above  stall).  The  nominal  potential  flow  for  an 
unstalled  airfoil  was  used  in  all  three  cases.  Results  obtained 
are  shown  in  Figure  5,  where  the  variations  between  separation 
and  reattachment  of  Ue,  Si  S * the  ordinate  where  u = 0 are 
compared . 

It  should  be  noted,  first,  that  converged  solutions  were  obtained 
for  all  three  cases  without  difficulty.  The  case  with  Cy  equal 
to  1.57  does  exhibit  a tendency  for  failure  to  reattach,  in 
that  near  reattachment  the  slope  of  the  ordinate  where  u = 0 is 
about  half  of  what  was  obtained  with  lower  values  of  . Never- 
theless, a converged  solution  was  obtained  with  a loading  about 
11  percent  greater  than  that  required  to  cause  bubble  burst, 
strongly  suggesting  that  failure  to  reattach  is  not  the  burst 
mechanism. 

Next,  an  attempt  was  made  to  determine  directly  whether  bubble 
burst  is  due  to  reseparation  of  the  turbulent  boundary  layer 
just  downstream  of  reattachment,  where  strong  interaction  between 
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Figure  4 - Variation  of  maximvun  lift  coefficient  with  Reynolds 
number  for  the  NACA  0012  airfoil  section. 


(X  -X  g)/c 

Figure  5 - Comparison  of  results  of  bubble  analyses  for 

three  different  lift  coefficients,  0012  airfoil 
Re  = 2 X 10°. 


the  viscous  and  inviscid  flows  must  still  be  taking  place. 

The  formulations  employed  to  analyze  the  bubble  shear  layers 
could  not  be  used  directly  to  analyze  the  turbulent  boundary 
layer,  because  the  turbulent  wall  shear  is  not  properly  taken 
into  account  in  those  relations. 

To  represent  the  turbulent  boundary  layer,  it  was  first  decided 
to  employ  the  finite-difference  method  of  Reference  5 and 
perform  a direct  iteration  between  viscous  and  inviscid  flows 
downstream  of  reattachment,  while  retaining  the  indirect  itera- 
tion procedure  for  the  bubble.  The  finite-difference  turbulent 
boundary  layer  code  was  successfully  incorporated  in  the  bubble 
analysis  program.  A displacement  thickness  variation  down- 
stream of  reattachment  was  then  obtained  from  an  initial 
estimate  of  the  Ug  variation  there,  for  an  airfoil  loading 
somewhat  below  stall.  Attempts  to  iterate  were  then  made, 
but  unfortunately  the  solutions  were  strongly  divergent,  the 
difficulty  apparently  stemming  from  extreme  sensitivity  of 
the  solution  to  the  pressure  gradient  just  downstream  of 
reattachment . 

It  was  then  decided  to  abandon  the  finite-difference  approach 
and  instead  use  an  integral  formulation  and  an  iterative 
procedure  analogous  to  the  one  used  for  the  bubble  itself.  The 
integral  relations  developed  for  the  turbulent  boundary  layer 
are  outlined  in  the  Appendix.  These  relations  were  success- 
fully incorporated  in  the  bubble  analysis  code,  and  again  a 
solution  was  obtained  for  the  boundary  layer  downstream  of 
reattachment  for  an  initial  assumed  pressure  distribution.  An 
iterative  procedure  was  then  set  up  which  is  exactly  analogous 
to  the  one  used  for  the  bubble.  That  is,  the  differential 
equation  for  g*  * was  used  to  generate  a linear  integro- 
differential  equation  for  the  flow  perturbation.  This  equation 
was  solved  in  exactly  the  same  manner  as  for  the  bubble,  again 
using  a trigonometric  series  to  represent  the  flow  perturbation. 
Unfortunately,  this  approach  was  not  successful  either,  the 
solution  for  the  flow  after  the  first  iteration  being  completely 
unrealistic.  The  cause  of  the  difficulty  is  not  presently 
clear.  A simple  error  in  the  formulations  or  coding  cannot 
be  ruled  out,  but  the  equations  and  program  were  carefully  checked. 
Another  possibility  is  the  method  used  to  integrate  through  the 
singularity  just  downstream  of  reattachment.  The  same  procedures 
as  were  used  in  the  bubble  analysis  were  employed  for  that 
singularity.  The  solution  appears  to  depend  strongly  on  the 
gradient  in  6"  * at  the  singularity  obtained  by  these  procedures. 

It  may  be  that  the  boundary- layer  equations  are  sufficiently 
different  from  those  of  the  shear  layer  that  a revised  approach 
is  required  to  integrate  through  the  singularity. 
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It  does  appear,  since  the  bubble  analysis  was  successful,  that 
a solution  for  the  turbulent  boundary  layer  in  the  vicinity  of 
reattachment  can  be  obtained  using  an  integral  formulation  and 
an  indirect  iteration  procedure.  Unfortunately,  the  limitations 
in  the  scope  of  this  study  precluded  further  pursuit  of  that 
approach . 

It  can  be  tentatively  concluded,  in  any  case,  that  the  specific 
mechanism  for  bubble  burst  is  not  failure  to  reattach,  but 
rather  must  involve  the  downstream  turbulent  boundary  layer. 
Further  analysis  should  be  directed  first  to  accurately  defining 
the  flow  in  the  boundary  layer  just  downstream  of  reattachment, 
taking  strong  interaction  effects  into  account. 
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APPENDIX 


Relations  for  Analyzing  the  Turbulent  Boundary  Layer 


The  procedure  for  analyzing  the  turbulent  boundary  layer  was  developed 
to  provide  an  analytic  form  which  is  compatible  with  the  relations 
used  for  the  shear  layer.  Specifically,  a pair  of  ordinary  first- 
order  differential  equations,  with  g*  and  a profile  parameter  as 
dependent  variables,  was  derived. 

Starting  from  the  momentum  integral  and  first  moment  equations, 
it  was  decided  to  employ  the  same  basic  approach  as  that  of  Albers 
(Reference  16)  whereby  the  equations  are  partially  decoupled  from 
the  local  pressure  gradient,  except  for  two  departures  from 
Albers'  approach.  Specifically,  rather  than  evaluate  dH/dx 
and  dJ/dx  separately,  it  is  assumed  that  dj/dx  = (dJ/dH) (dH/dx) , 
where  dJ/dH  is  a function  of  profile  parameter  a and  the  friction 
coefficient  Cf.  Also,  a relation  for  Cf  was  adopted  which  closely 
approximates  that  of  the  Ludwieg-Ti liman  correlation  (Reference 
16)  near  separation,  namely 


Cf  = .041  a Rg 


-.268 


where  the  profile  parameter  a is  defined  by 

a = H _i  = -H_ 

Hg  .429  ^ 

and  RgQ  is  local  momentum  - thickness  Reynolds  number. 
The  specific  relations  derived  are  then  as  follows. 


d ^ * 
dx 

da  _ 
dx 


where 


= (HdJ/dH-J)  ( [3J-(2H+l)dd/dH]  ^du,/dx'j 

^ ( HD-f2j+J(l-H)  dUp/dx  \ 

J-*  (HdJ/dH-J)  [ Ug  J 

f2  = Cf/2 
H = .429(a+l) 

J = 3H-I+I3/I2 

dJ/dH  = 3 [l-(1.7f+1.2994/^)/K  ] 
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